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Accurate quartet- and doublet-state potential-energy surfaces for the interaction of a hydrogen atom and an 
oxygen molecule in their ground states have been determined from an ab initio calculation using large-basis 
sets and the internally contracted multireference configuration interaction method. These potential surfaces 
have been used to calculate the H 0 2 electron-spin-exchange cross section; the square root of the cross section 
(in a 0 ), not taking into account inelastic effects, can be obtained approximately from the expressions 
2.390£” 1/6 and 5.266-0.708 log l0 (£) at low and high collision energies E (in E h ), respectively. These func- 
tional forms, as well as the oscillatory structure of the cross section found at high energies, are expected from 
the nature of the interaction energy. The mean cross section (the cross section averaged over a Maxwellian 
velocity distribution) agrees reasonably well with the results of measurements. ^\y / 1 S 
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I. INTRODUCTION 

Spin-exchange collisions are of interest for studies of 
gases of the upper atmosphere [1] and interstellar space [2]. 
Accurate interaction energies [3,4] have allowed accurate 
calculations [5-7] of the spin-flip cross sections for the col- 
lision of hydrogen atoms; the theoretical results agree well 
with measured data [8]. On the other hand, the comparison 
for H-0 2 is considerably less satisfactory; spin-flip-scattering 
calculations with accurate interaction energies have been rec- 
ommended [9,10] to help resolve differences in measured 
data. 

Walch and co-workers [11 — 13] have reported extensive 
ab initio calculations to define the ground-state H0 2 
potential-energy surface; references to earlier work are con- 
tained therein. Unfortunately, their results do not include the 
potential data for large H-0 2 separation distances that are 
required to determine spin-flip cross sections, especially at 
lower collision energies. We have calculated accurate H-0 2 
potential-energy surfaces for the doublet and quartet states, 
corresponding to an interaction of the atom and molecule in 
their ground states, for a broad range of separation distances. 
We have extended these results with the proper long-range 
forces to determine transport collision integrals [14(a)] and 
spin-flip cross sections for a broad range of energies. A de- 
scription of the molecular- structure calculations and the con- 
struction of the potential-energy surfaces is presented along 
with results in Sec. II. Approximations to the spin-flip cross 
sections and an outline of the scattering calculation can be 
found in Sec. Ill; the spin-flip cross sections from the scat- 
tering calculation are presented and compared with approxi- 
mations and measured data in Sec. IV. Concluding remarks 
are contained in Sec. V. 

II. DETERMINATION OF THE POTENTIAL-ENERGY 
SURFACES 

The potential-energy curves are determined from ab initio 
calculations using the complete-active-space self-consistent- 


field (CASSCF)-intemally contracted multireference con- 
figuration interaction (ICMRCI) method [15,16]. The calcu- 
lations are performed in C s symmetry with the oxygen 2 p 
orbitals and the hydrogen Is orbital active. The oxygen 2s 
orbitals are correlated in the Cl calculation, but are con- 
strained to be doubly occupied in all reference configura- 
tions. A multireference analog of the Davidson correction 
[17,18] (+Q) is used t0 estimate the effect of higher excita- 
tions. 

The one-particle basis sets employed are the augmented 
correlation-consistent polarized- valence basis sets by Dun- 
ning and co-workers [19,20], The triple-^ (TZ) basis set was 
employed for extensive calculations that could define a com- 
plete and nearly complete potential-energy surface for the 
doublet and quartet states, respectively. A smaller number of 
calibration calculations were performed employing the larger 
quadruple-^ (QZ) basis set. Basis-set superposition errors 
(BSSE’s) were determined using the counterpoise method 
[21]. The ab initio energies corresponding to the QZ calcu- 
lation will be tabulated elsewhere [14(b)], both the TZ and 
QZ results are, however, available from the authors. 

The coordinate system for the present work is specified by 
the separation distance r of the H atom from the center of 
mass of 0 2 and an angle y between a line from H to the 
center of mass of 0 2 and the 0 2 symmetry axis passing 
through the nuclear centers. Making a rigid-rotor approxima- 
tion for 0 2 , the 0-0 separation distance tqq is fixed at the 
equilibrium value 2.28a 0 . For scattering calculations, a bet- 
ter choice for roo would be the expectation value for the first 
vibrational state; however, since the size of the spin-flip 
cross section is primarily determined by the interaction en- 
ergies at large r (see the analysis of Sec. Ill below), one 
expects that the correction to the spin-flip cross section from 
this slightly larger value for Tqo will not be significant, i.e., 
within the uncertainty of the present scattering calculation. 
The potential energy curves V(r, y) derived from the present 
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FIG. 1 . H-0 2 potential-energy curves for the doublet state for 
various values of the angle y. All curves were constructed from 
spline fits to the discrete data. 

calculation are shown in Figs. 1 and 2, for the doublet and 
quartet states, respectively. 

The doublet-state potential-energy curves for six values of 
y (0°, 30°, y 0 , 60°, 75°, and 90°) were determined from QZ 
calculations. The angle y 0 ^54.1° is defined by the condition 

P 2(cosy 0 ) — 0) 

where P 2 is a Legendre polynomial of order 2; the potential 
energy at this angle is of interest for scattering calculations 
(see Sec. IV below) since it becomes a good approximation 
to the spherically averaged potential energy at large r. These 
curves and the results of the TZ calculations were then used 
to construct curves for additional values of y. Corrections to 
the TZ results were obtained from least-squares fits to the 
energy differences A V(r,y) y between the energies from the 
QZ and TZ calculations, using Legendre polynomials 
P n (c os y) with n restricted to even values by symmetry. For 
large r, where the correction is relatively large because of 
BSSE’s for the TZ results, we found that the corresponding 
fit to the logarithm of AV(r,cos y) yielded higher accuracy 
for the predicted corrections. Some of the data points at 
about 3 .5a 0 for the potential-energy barriers, shown in Fig. 1, 
for small y had to be determined directly from QZ calcula- 
tions; here the fitting is not accurate because of the large 
variation in AV as y approaches zero. In addition, the TZ 
results were used to construct the potential-energy walls at 
small r; here the improvement from the QZ calculation is 
expected to have only a negligible effect on the results of the 
present scattering calculation. Furthermore, taking into ac- 
count the rigid-rotor approximation for 0 2 of the present 
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FIG. 2. H-0 2 potential-energy curves for the quartet state for 
various values of the angle y. All curves were constructed from 
spline fits to the discrete data. 

calculation, one expects that these QZ calculations will most 
likely not yield a meaningful change in the scattering results. 

The quartet-state potential-energy curves were con- 
structed from the results of the QZ calculation for the same 
angles used for the doublet case, excluding 60°. The curves 
for y= 15° and 45° were obtained by the fitting procedure for 
the energy corrections described above for the doublet case. 
Because of the more uniform behavior of the quartet 
potential-energy surface (in contrast to the complex nature of 
the corresponding doublet-state curves arising from avoided 
curve crossings [11-13]) fewer curves at different values of 
y are required to define a complete potential energy surface. 
The curves for remaining values of y shown in Fig. 2 were 
determined from least-squares fits to the potential-energy 
data for the above seven values of y, using the P„ and the 
procedures described in the preceding analysis. As in the 
doublet case, the repulsive walls for the quartet state were 
also constructed from the results of the TZ calculation at 
small r. 

The coefficients of the various polynomial fits that are 
described above, are tabulated in Ref. [14(b)] to facilitate the 
determination of V(r, y) for values of y not covered in this 
work. The accuracy of the potential data from the fitting 
procedures was determined by comparison with the results of 
QZ calculations. At large r, for example, we found that for 
the quartet state the predicted values of V(r,75°) determined 
from fits to In A V using only the data for the other four 
angles (0°, 30°, y 0 , and 90°) agree with the corresponding 
energies from the QZ calculation to within 1 % in the van der 
Waals region (r^5a 0 ). Another confirmation of accuracy 
will be pointed out in the following paragraph. 
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FIG. 3. Energy difference between the potential energies of the 
quartet and doublet states for selected values of the angle y. A pair 
of curves, one that represents the QZ (solid line) and another that 
represents the TZ (dashed line) calculations, is shown for each 
value of y. 

At higher energies, the size of the spin-flip cross section is 
(to a first approximation, see Sec. Ill below) determined by 
the magnitude of the difference between the potential ener- 
gies of the doublet and quartet states at large r. The energy 
difference obtained from the QZ calculation is compared 
with the corresponding results for the TZ calculation in Fig. 
3. The close agreement of the two sets of curves suggests 
that the present calculation yields a nearly converged (accu- 
rate) value for this quantity and therefore that the present 
potential curves should yield accurate scattering cross sec- 
tions. Furthermore, the agreement shown by the curves for 
15° and 45° provides additional support for the accuracy of 
the fitting procedure described above. 

The ab initio curves have been extended at large r with 
the long-range expansion 

V(r,y)-> X % [1 +r 2n / > 2 (cosy)] (2) 

n = 3 r 

where C 2n is the isotropic dispersion coefficient and the pa- 
rameter r 2 „ specifies relative anisotropy. For the leading 
term, we take the value 19.16£^ao for C 6 determined by 
Zeiss and Meath [22] from oscillator strength, photoabsorp- 
tion, and scattering data and the value 0.21 1 for T 6 obtained 
by Langhoff, Gordon, and Karplus [23] from optical disper- 
sion absorption and scattering data. Noting that the ratios of 
the dispersion coefficients for interactions involving hydro- 
gen are about the same (e.g., see Ref. [24]), we can obtain 
estimates of the two higher-order dispersion coefficients us- 


ing the values of C 2n calculated by Meyer (as reported in 
Ref. [25]) for H-H 2 ; i.e., 8.813E*a$, 162.3 E h al, and 
3999E h a' 0 °, for n =3-5, respectively. Hence, finding that the 
ratio of C 6 for H-0 2 compared to H-H 2 is 2. 17, we obtain the 
estimates 352 E h a% and 8678 E h a' 0 ° for n= 4 and 5, respec- 
tively, for the present work. 

Values for the higher-order (n=6-8) coefficients were 
obtained from recursion relations [26,27] 

^2n+4 = (^2n + 2 ^2n) ^2n-2 • (3) 

We found that the long-range interaction energy calculated 
from the above values of C 2n and Eq. (2) for y 0 [where the 
second term within the square brackets of Eq. (2) vanishes] 
agrees well with the corresponding results from the QZ cal- 
culation for this angle at large r. Furthermore, comparing the 
long-range interaction energy with the results of the QZ cal- 
culation for other angles, we found that including the aniso- 
tropic contribution for n— 4 with the value 1.15 for 1 8 im- 
proved the agreement of the values of V(r,y) from Eq. (2) 
with the corresponding ab initio results. 

HI. THE SCATTERING CALCULATION 

A general formulation of the scattering for collision part- 
ners, such as the hydrogen atom and the 2-state oxygen mol- 
ecule of this study, has been developed by Launay [28] using 
formalism for the body-fixed system. We have adopted the 
sudden approximations of Parker and Pack [29] to calculate 
the scattering in the body-fixed system; the spin-flip cross 
sections from this approach should be sufficiently accurate, 
provided the collision energy E is not too low, Applying both 
the energy and centrifugal sudden approximations of Ref. 
[29], we find that the the spin-flip cross section (for homo- 
nuclear molecules) can be obtained from 

f tt/2 

tr SF (£)= o- SF (£,y)siny dy, (4) 

Jo 

where the cross sections cr SF (£\y) are determined for an ori- 
entation specified by the angle y that remains fixed during 
the collision. Thus a sf (E y y) can be calculated from a central- 
field formulation of the scattering [9], i.e., from 

X 

o- SF (£,y)= 7 T 2 (21+ \)sm 2 [tf(E,y)- itf(E,y)], 
k /=o 

(5) 

where / is the angular momentum quantum number and k is 
the wave number. The scattering phase shifts rf^{E y y) are 
calculated from the interaction energies for the quartet and 
doublet states (/jl ~ 4 and 2, respectively) at each fixed value 
of y. 

At this point, we shall pause to examine the behavior of 
the spin-flip cross section, inferred from the interaction en- 
ergies described in the preceding section and the above scat- 
tering approximations. 

We can obtain a low-energy estimate of cr SF by generaliz- 
ing the analysis of Rapp and Francis [30] for charge- 
exchange collisions, which is based on the results of Gio- 
mousis and Stevenson [31] for a long-range polarization 
force. Following their arguments for scattering in an attrac- 
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tive long-range force field, only the penetrating collisions 
with impact parameters b^b oy the impact parameter for 
classical orbiting, can contribute to spin exchange. Taking 
the probability of spin exchange to be 1/2 for these collisions 
[i.e., the average value of sin J in Eq. (5) for or SF ], using the 
semiclassical relation b = (l + 1/2)//:, and replacing the sum- 
mation of Eq. (5) by an integration, we obtain 


b db= jvbl- (6) 

Jo 

For a given value of £, b Q specifies that £ is equal to the 
value of the maximum of the potential-energy barrier exhib- 
ited by the effective potential energy V e (r,b)> i.e., b 0 satis- 
fies 

b 2 0 

E-V e (r,b 0 ) = E-E-^-V(r) = 0, (7) 


x v e (r,b 0 ) = 0. (8) 



Taking the long-range form 

V(r)=-^ (9) 

and combining Eqs. (7) and (8), we obtain 

pip \ i 

^ = 2\2 _1 J ( C ^ /£ ) 2/V - 00) 

The spin-flip cross section can be readily obtained by com- 
bining relations (6) and (10); for the case i/=6, we find 

tt(C 6 IE ) ]I 3 . (11) 

The nearly linear behavior and common slope of the en- 
ergy difference curves shown in Fig. 3 for the values of y 
corresponding to the leading contributors to the scattering 
indicates that an analytical approximation can also be con- 
structed for the spin-flip cross section at high energies. Ap- 
proximations to the cross sections of the form (5) have been 
developed [32,33] for an energy difference that can be rep- 
resented by an exponential at large r; their application to the 
present potential-energy data is illustrated in the following 
section. For the purposes of the present paper, however, we 
point out that the major contributions to spin flip comes from 
orientations where the corresponding cross sections can be 
represented [33] by a linear expansion in In (£); i.e.. 


FIG. 4. H-0 2 spin-exchange cross section. The dotted and 
dashed lines represent the least-squares fits to the discrete data; the 
fits are based on the functional forms of Eqs. (11) and (12), respec- 
tively. 

[o- SF {£ (r )] ,/2 =A(y) + fi( r )ln(£). (12) 

Provided the variation from the second term on the right- 
hand side of Eq. (12) is relatively small, it follows from Eq. 
(4) that [<j sf (£)] 1/2 can also be represented by a linear ex- 
pansion in ln(£). 

At low energies, the phase shifts for the present work 
were obtained by direct numerical solution [34] of the Schro- 
dinger equation. At higher energies they were obtained from 
a semiclassical method that includes a uniform approxima- 
tion [35,36] that accounts for the quantum-mechanical ef- 
fects arising from a potential-energy barrier. The integration 
for Eq. (4) was accomplished with Gaussian quadrature. The 
results of the scattering calculation are presented and dis- 
cussed in the following section. 

IV. SPIN-FLIP CROSS SECTIONS 

The spin-flip cross section from the present calculation is 
displayed in Fig. 4 along with analytical least-squares fits 
based on the function forms of the preceding section. The 
energy dependence of tr SF from Eq. (11) fits the low-energy 
data well. Taking the value for C 6 for H-0 2 selected in Sec. 
II, we find that Eq. (11) yields 

(7^~2.8£“ 1/6 , (13) 

where E is in units of E h and cr SF is in Comparing this 
result with the fit shown in Fig. 4, we find that Osp is over- 
estimated by about 17%. 
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As mentioned above, the potential-energy curve for y 0 is a 
good approximation to the spherically averaged potential en- 
ergy V(r) at large r. In previous work [37,38] we have found 
in some cases that elastic-scattering cross sections deter- 
mined from V(r) can be a good approximation to the corre- 
sponding result obtained from the complete potential-energy 
surface. From a least- squares fit to the present QZ data (see 
Fig. 3), we obtain 

V A (r,y 0 )-V 2 (r,y Q )=n3\ exp( - 1 .7377r); (14) 

taking this result (14) and applying the approximate methods 
[32,33] pointed out in the preceding section, we find 

[<r SF (£)] 1/2 ~ 5.810-0.882 lo gl 0 (£), (15) 

where the units are the same as specified for Eq. (13) above. 
This approximation yields a line that lies slightly below 
(about 2% lower than) the corresponding calculated data on a 
semilogarithmic plot such as shown in Fig. 4; i.e., the term 
corresponding to A(y) of Eq. (12) would require only an 
increase of about 0.1 a 0 to reach agreement with the results 
from the scattering calculation. Comparing the values from 
Eq. (15) with those from the high-energy least-squares fit to 
the calculated data shown in Fig. 4, we find that cr^p for y 0 is 
about 12% higher than the mean representing the calculation 
for the complete potential-energy surface. 

Note the oscillatory structure exhibited by the results of 
the scattering calculation shown in Fig. 4 at higher energies. 
These oscillations are expected from the behavior of the po- 
tential curves shown in Figs. 1 and 2; the presence of poten- 
tial barriers produces extrema in the difference V 4 — V 2 . Cal- 
culations [36,39] have shown that these extrema cause 
oscillations in the cross section. In fact, we find such oscil- 
lations in (r SF (£,y) for those values of y for which the cor- 
responding potential curves have barriers. 

The mean spin-flip cross section, obtained by averaging 
over a Maxwellian velocity distribution, i.e., 

& sf (T) = (kT)- 2 f a SF (E)E e\p(-E/KT)dE, (16) 
J 0 

where k is the Boltzmann constant and T is the kinetic tem- 
perature, is of interest for applications. The values of cr SF (T) 
from the results of the present calculation are shown as a 
function of T in Fig. 5 along with an analytical approxima- 
tion to facilitate future applications of these results. 

We also compare the present results with measured data 
[40-42] using the corrections of Turner, Snider, and Fleming 
[10] in Fig. 5. The value from the most recent room- 
temperature measurements by Anderle et al. [42] agrees well 
with the result from our work. The lowest-temperature mea- 
surement of Gordon et al. [41] also agrees well with present 
results; on the other hand, the higher-temperature data fall 
off more rapidly [more like the T 1/3 low-temperature ex- 
trapolation expected from <t sf (£') with the energy depen- 
dence of Eq. (11); see Fig. 5]. 

V. CONCLUDING REMARKS 

We have determined the rigid-rotor H-0 2 potential-energy 
surfaces that cover the complete range of separation dis- 
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FIG. 5. H-0 2 mean spin-flip cross sections in A 2 from the re- 
sults of the present calculation (solid line) and from measured data 
(the vertical lines represent the error bars for the corresponding data 
point). The dashed line represents an analytical approximation to 
the present results. The dotted line was determined from the value 
(21.85 A 2 ) of spin-flip cross section at 100 K from the present 
results and the temperature dependence of cr SF ( T) that is obtained 
from Eq. (16), if tr SF (E) is proportional to £ -1/3 as in Eq. (11). 

tances including the long-range interaction energies required 
for calculating scattering cross sections. Our analysis indi- 
cates that the energy difference between the quartet and the 
doublet states obtained from our results is very accurate at 
large separation distances and consequently that the present 
potential-energy data should allow an accurate determination 
of the spin-flip cross section. The present scattering results 
should be accurate at the higher energies, but at lower ener- 
gies the present potential data merit a more accurate treat- 
ment than that of the present calculation. For example, at 
large impact parameters, one might follow the approach of 
Stallcop [43] that retains only the energy sudden approxima- 
tion. Another approach has been mentioned in the beginning 
of Sec. Ill above. 

We have compared low- and high-energy approximations 
for as F with the scattering results of the present work; this 
comparison should contribute to the understanding of the 
physics of collision-induced spin-flip processes and the ap- 
plication of the approximations for estimating <t sf (E) from 
limited potential -energy data. 
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